Tunneling into clean Heavy Fermion Compounds: Origin of the Fano Lineshape 
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Recently observed tunneling spectra on clean heavy fermion compounds show a lattice periodic 
Fano lineshape similar to what is observed in the case of tunneling to a Kondo ion adsorbed at the 
surface. We show that the translation symmetry of a clean surface in the case of weakly correlated 
metals leads to a tunneling spectrum given by the superposition of the local weighted density of 
states of all energy bands involved, which does not have a Fano lineshape. In particular the spectrum 
will show any hybridization gap present in the band structure. By contrast, in a strongly correlated 
heavy fermion metal the heavy quasiparticle states will be broadened by interaction effects. The 
broadening grows as one moves away from the Fermi surface, up to a value of the order of Tk , the 
Kondo scale. We show that the hybridization gap is completely filled in this way, and an ideal Fano 
lineshape of width Tk results, similar to the impurity case. We also discuss the possible influence 
of the tunneling tip on the surface, in (i) leading to additional broadening of the Fano line, and (ii) 
enhancing the hybridization locally, hence adding to the impurity type behavior. The latter effects 
depend on the tip-surface distance. 



Introduction. - Recent progress in scanning tunneling 
spectroscopy (STS) techniques has made it possible for 
the first time to measure the differential conductance, 
dl/dV, in heavy-fermion compounds [l], 0]. These ma- 
terials arc governed by the strong Coulomb interaction 
on the f-electron orbitals, hybridizing with a conduc- 
tion band. Provided the f-energy is sufficiently below 
the Fermi energy and the hybridization strength is weak 
enough, magnetic moments will appear at the f-sitcs. The 
Kondo effect will screen the local moments at a tem- 
perature below the Kondo temperature Tk @, 0] • The 
ensuing Fermi liquid state is characterized by large effec- 
tive masse ratios m*/m . Out of this heavy Fermi liq- 
uid a variety of interesting phenomena may emerge, such 
as magnetic order, quantum phase transitions and the 
associated non-Fermi liquid behavior, as well as uncon- 
ventional superconductivity [H[ . Their microscopic origin 
likely lies in the competition between the exchange inter- 
action of the magnetic moments and their screening by 
conduction electrons Iff, though no theoretical consensus 
has emerged as yet [2J. STS experiments, by providing 
insight into the local electronic structure [l|, [2| of heavy- 
fermion materials, might hold the key to understanding 
their complex properties, by providing complementary 
information. The theoretical challenge in the interpre- 
tation of the differential conductance in Kondo lattice 
systems [H, Q lies in the proper treatment of the strong 
correlation effects governing these systems. Whereas tun- 
neling into single Kondo impurities IB" 12 1 is relatively 
well understood, there is a serious problem of interpre- 
tation of tunneling into states on a clean crystal surface, 
which has the lattice periodicity. In the single impurity 
case the quantum interference between electrons tunnel- 
ing from the STS tip into the conduction band and into 
the magnetic f-electron states is essential. As a conse- 



quence, the lineshape of the dl/dV spectrum shows an 



asymmetric form, as first discussed by Fano [18(. It is 
rather surprising that a similar spectrum is found, peri- 
odically continued along the surface, in the clean case. 
Within the framework of electron band theory one would 
rather expect the spectrum to show the superposition of 
local density of states contributions from the different 
electronic bands involved. Recent attempts to calculate 
the tunneling spectrum by treating the problem of strong 
correlations within pseudoparticle mean field theory bear 
out this expectation [H, Q. The resulting spectra are 
characterized by two peaks belonging to the two heavy 
bands, separated by a hybridization gap. As we will show 
below, a mapping of the problem onto a non-interacting 
system, as done in 1/N theories is not sufficient to cap- 
ture the dominant many-body effects at energy scales of 
the order of Tk ■ 

In this paper, we address this issue within Fermi liquid 
theory extended to somewhat higher energy scales. It is 
convenient to start with the lattice Anderson model in 
the low energy domain where heavy quasiparticles (QP) 
are well defined. Our main point is that the QP width 
r at energy E, which varies as T oc E 2 /Tk in the Fermi 
liquid regime (we use units with H = krj = f), will grow 
up to r oc Tk at the scale E = Tk, which is larger than 
the hybridization gap, given by A cx (V/e° F )TK « Tk, 
where V is the hybridization amplitude and e F is the 
bare Fermi energy of the unhybridized conduction band. 
Consequently, the local f-electron density of states (DOS) 
takes the form of a Lorentzian of width Tk, just like 
in the impurity case, and the Fano lineshape arises as 
usual. We will first consider the case of a weakly corre- 
lated metal, to make our point that in this case no Fano 
lineshape is expected. 

While the translation invariance is preserved in an 
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ideal noninvasive tunneling experiment, in reality the 
translation invariance may be broken by the local influ- 
ence the tip may exert on the surface. Such an effect has 
been recently invoked to explain the somewhat distorted 
energy spectrum observed in tunneling spectra of a topo- 
logical insulator [f9||. We discuss two ways by which the 
presence of the tip may change the tunneling spectra. 
Both effects will depend on the tip-surface distance. 
Tunneling into weakly correlated metals. - We consider 
tunneling into a metal characterized by two hybridized 
bands, as described by the Hamiltonian 
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Hhyb 



Ho + Hhyb + H t 
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where cjL., f^,pt create an electron of spin a in a Bloch 
state of momentum k, a localized f-orbital at site i or the 
level at the tip of the tunneling electrode, respectively. 
The operator n/j CT = ftaha counts the number of elec- 
trons with spin a on the local f- level, and e/ is the energy 
of the f-level (which is position and momentum indepen- 
dent). In Eq. ([1]) we only take into account tunneling 
into the orbitals directly under the tip at Ro = 0, which 
is the only source for breaking of the translational in- 
variance, otherwise preserved in the Hybridized system. 
Assuming the tunneling electrode and the metal to be in 
thermal equilibrium (in the limit of vanishing tunneling 
current), their chemical potentials differing by eV, where 
V is the applied bias voltage, the tunneling current to 
lowest order in the tunneling amplitudes t Cl tf (taken to 
be real valued) is given by 



I(V) = ^[<kjN t (u-eV)[f(u-eV)-f(u)]x 

x ${t 2 c G cc (u) +t}G ff (oj) + 2t c t f G cf (w)} (2) 

where G a t,(u>) are the advanced local single particle 
Green's functions of H wc = Ho + Hhyb and N t (uj) is 
the STM tip DOS. The Hamiltonian H wc may easily be 
diagonalized to yield two hybridized bands with Bloch 
energies 



and density of states 



N +t -(u>) = / d J H(w 
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FIG. 1: Density of states for a two-band model of Eq.[T]with- 
out the QP width (dashed line) and with the QP width of 
Eq. ll2l (solid line) with Tk — 130K (see text for other parame- 
ters). The inset shows the DOS with QP width (circles) and a 
fit to the Fano lineshape (solid line), yielding a Fano linewidth 
of ~ 20mcV, in good agreement with experiment P, @1- For a 
qualitative comparison, in the inset we plot the experimental 
Fano lineshape, taken with permission from Ref. [l| . 



The density of states exhibits a hybridization gap at lj = 
ef of width A = 2V . The Green's functions are found 
as 



G a b(uj) 



E 

v=+,- 



ih,k 



2Ae k 



V 



2Ae k 



2Ae k 



(5) 



where Ae k = — ejT. While the coherence fac- 
tors a v ab k shift the weight in the partial spectral func- 
tions 3{G a f,(w)} somewhat, compared to the band den- 
sity of states N+,-(u), the hybridization gap will re- 
main. In Fig. [T] we show the differential conductance 
dl/dV calculated for a three-dimensional parabolic band 
of bandwidth 2eV (crossing the Fermi energy at \k\ = 
0.8667r/ao), Hybridization strength V = O.leV, the tun- 
neling amplitude ratio tj/t c = 0.15 and f-level energy 
e f = ~ lOeV without any finite level width (dashed line). 
The shape of the curve is similar to what has been found 
in Refs. [1, Q within pseudoparticle mean field theory of 
the Kondo lattice model. These results are nowhere near 
what is observed in experiments on strongly correlated 
heavy fermion compounds. We stress again the main con- 
clusion from this section: in the absence of finite f-level 
lifetime and breaking of translation invariance, the DOS 
cannot exhibit a Fano lineshape. Therefore, understand- 
ing the origin of either finite QP lifetime or a breaking 
of translational invariance of the system by the influence 
of the tip on the surface is crucial to understanding the 
Fano lattice observed in recent experiments. 
Tunneling into heavy fermion compounds. - A micro- 
scopic Hamiltonian believed to capture the physics of 
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heavy fermion compounds is the lattice Anderson model 



Hal = H u 



U J2 



n M n f*l 



(0) 



where U is the Coulomb interaction matrix element, 
assumed to be much larger than the conduction band 
width. 

As above the single particle Green's functions as well as 
the self-energy form 2x2 matrices in the conduction elec- 
tron (c) , f-electron (f) space. It is worth noting that the 
only self energy element is Titf{uS) . Expanding Ejy(ct>) 
about the Fermi energy u> = (neglecting its weak mo- 
mentum dependence), we end up with a QP description 
in terms of the QP weight factor z , the shifted position 
of the f-level e/ and the QP line width (the inverse of the 
QP life time) 7 , which are defined in terms of E//(w) 
by 



z- 1 = l-—StE f/ {w)\ , 
e f = z[e f + RE ff (0)] 
7 (£ k ) = *££//(£*), 



(7) 
(8) 
(9) 



The complex-valued QP energies of the two hybridized 
bands are defined by 



1 



£k = ^{e f + e u -i 1 ±^{e f -e k -i 1 ) 2 +4zV 2 



(10) 



The hybridization gap follows (neglecting the QP width) 
as A = zV 2 ( \ 1 + 1 — - — r) i where G is a vector at 

M e k=o |£k=G|' ' 

the edge of the Brillouin zone. In the heavy fermion 
regime, assuming that the Fermi energy (-Ek F = 0) in- 
tersects the lower band and observing >> 1 , and 
£k F >> v^l^kl one ma Y expand about the Fermi 
momentum to get E k = (k — kp)vp m »"g k ) , where 



we defined the effective mass ratio 



Here tk F is the value of the bare conduction band en- 
ergy at the Fermi surface, which is assumed to be close 
to the upper band edge , so that e^ F sa W/2, where W 
is the conduction band width. Approaching the heavy 
fermion regime from high temperatures the renormaliza- 
tion is observed to occur below a scale Tk , usually called 
the lattice Kondo temperature. We assume Tk to be ap- 
proximately equal to the Kondo temperature of a single 
Kondo ion. 

It follows from the Pauli principle that the QP width 
tends to zero in the limit temperature T and excitation 
energy E k — > , as given by 



T(E k ) = A[El + (ttT) 2 ] « £ k 



(11) 



As long as this condition is satisfied, one may employ 
the formalism of Fermi liquid theory. There are indica- 
tions from both , theory and experiment (for a recent 



determination of T from ESR data see 2l|) that the 
prefactor scales with the inverse Kondo temperature as 
A ss 1/Tk- At the energy scale Tk we therefore have 
T(Tk ) ~ Tk ■ Beyond that scale the lattice coherence is 
suppressed by inelastic processes and the metal behaves 
like a crystal of independent Kondo ions. The coher- 
ence scale of the single Kondo ion is set by the local 
spin relaxation time, which is known to be T s rj Tk at 
energy /temperature less or equal to Tk ■ At higher en- 
ergy the relaxation rate is given by T S (E) sa g 2 (E)\E\, 
where g{E) « 1/ \n(\E\/T K ), E > 2T K ■ At tempera- 
tures T << Tk one may therefore approximate the QP 
width by 



T(E k ) = El/T K , £ k < T K 
T(E k ) = \E k \/[l + H\E\/T K )} 2 



E k >T K .(12) 



The solid points in Fig. [T] denote the DOS calculated 
in the presence of the finite QP width of Eq. [T2l taking 
T K = 130K 0. The solid black line is a fit to the Fano 
lineshape , given by [HI p F ano = where e = ^Si. 

The fit yields q = 1.25 and T = 20.7 meV, in good agree- 
ment with the values obtained in experiment [H, Q . For 
a qualitative comparison, in the inset we plot the exper- 
imental Fano lineshape of Ref. [l[ . From calculations of 
the DOS and fitting to a Fano lineshape with various pa- 
rameters, we find the relation T ~ 2Tk to hold for a wide 
range of parameter values. 

In the experiment [lj the Fano width T also varies with 
the STM tip position, although the modulation is smaller 
than that of the Fano asymmetry parameter q. While 
for the latter the origin of the modulation is clear (since 
the tunneling amplitudes t c and tf depend on the tip 
position), a modulation of T is somewhat surprising. A 
possible interpretation is provided in the following. 
Tip-induced effects. - We point out two ways in which 
the tunneling tip may influence the properties of the sur- 
face and consequently the tunneling lineshape. First, the 
f-electron self-energy acquires a contribution originating 
from tunneling into the tip, 9Sj/ jt i p (a;) oc \tf\ Nt(u)), 
resulting in a modulation of the QP width with tip posi- 
tion. The fact that this is only a part of the self-energy 
accounts for the relatively weak modulation in P. 

A second and potentially more important effect is re- 
lated to the local potential generated by the tip at the 
surface underneath [l9j . This potential acts like an ad- 
ditional local hybridization between f- and c-states, and 
will give rise to Kondo impurity type behavior. More 
specifically, in addition to the interaction-induced broad- 
ening, there will be a tip-induced local scattering term 
added to the Hamiltonian of Eq. [TJ of the form 



Hup — Vti p c r _ / r= o + h.c. 



= Vt 



tip 



(13) 
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FIG. 2: Density of states of the two-band model Eq. Q] in- 
cluding the QP width of Eq. [12] and the the effect of the 
tip-induced hybridization, for various values of the tip poten- 
tial Vup ■ The inset shows the variation of the Fano width 
with increasing Vu p 



where r = is the tip position and Vu p is the hybridiza- 
tion strength, which will depend on the tip position as 
well as on other parameters such as tip height above the 
surface and operating voltage. It is simple to show (20| 
that this kind of hybridization, which breaks translation 
invariance, directly leads to a Fano lineshape even in the 
absence of c-f hybridization. In the presence of the hy- 
bridization and the intrinsic broadening described above, 
the DOS can be evaluated by solving the Dyson equation 
in the c-f : 

G kk , =Gl5(k,k')+Gl-Vu P -{l-V Up -xy 1 -G kl , (14) 

where G° is the 2x2 Green's function matrix \ = J^k ^fe 
and Vup = Vup ( ^ ^ J. In Fig. [2] we plot the DOS calcu- 
lated with Vup = 0, 5, 15, 25, 35 meV (shown by arrows in 
the figure). The lineshape seems to be hardly affected, 
however, a fit to the Fano lineshape shows that while 
the Fano asymmetry parameter q hardly changes (within 
10%), the Fano width T, depicted in the inset as a func- 
tion of Vtip, changes by 150%, as seen in experiment. 
We therefore conclude that while the dominant cause for 
modulation of q is the modulation of the tunneling am- 
plitudes, the modulation in T may be due to modulation 
of the tip-induced local scattering. 

Conclusion. In this paper we offer an interpreta- 
tion of the recent observation of a lattice periodic pat- 
tern of Fano-shaped tunneling spectra of heavy fermion 
metals, which appears to be very similar to what is ob- 
served at Kondo impurity atoms on the surface. First 
we demonstrate for the case of weakly correlated hy- 
bridized bands that the tunneling spectra reflect the lo- 
cal density of states, which does not show a Fano line- 
shape in general. We propose the intrinsic correlation 



induced quasi-particle line broadening as the origin of 
the Fano lineshape observed in tunneling experiments on 
strongly correlated metals. We model the quasi-particle 
energy width phcnomenologically by using Fermi liquid 
theory and Kondo physics. In fact the broadening is re- 
lated to the Kondo temperature via Eq. and yields a 
Fano linewidth T ~ 2T k , in agreement with recent exper- 
iments. 

In addition, we address the modulation in the Fano 
width r observed experimentally. While the lattice peri- 
odic modulation of the Fano asymmetry parameter q is 
immediately understood as a consequence of the modu- 
lation of the tunneling amplitudes between tip and f- or 
c-orbitals, the modulation of the width T requires a differ- 
ent interpretation, possibly related to the influence of the 
tip on the properties of the metal beneath it. Specifically, 
we show that a tip-induced local hybridization may gen- 
erate such a modulation. A direct consequence would be 
that the Fano parameters q and T would depend on STM 
tip parameters (height, voltage etc.). A dependence of q 
on tip parameters is rather expected, since is depends on 
the tunneling amplitudes t c and tf, which are unlikely 
to change in the same way with changing tip parame- 
ters, hence tf/t c should depend on the tip parameters. 
However, a dependence of T on tip parameters would be 
surprising if T is a purely intrinsic quantity, and would 
indicate that indeed the STM tip effect in non negligible. 
These predictions can be directly tested within current 
experimental setups. 

We acknowledge stimulating discussions with C. Sea- 
mus Davis, M. Graf, J. Fransson and P. Coleman. PW 
acknowledges partial support by the Theory Division of 
Los Alamos Laboratory, the Aspen Center for Physics 
and the DFG research unit " Quantum phase transitions" . 
AVB and YD acknowledge support by US BES and 
UCOP-TR-01. 



[1] A. Schmidt et aJ.„ Nature 465, 570 (2010). 

[2] P. Aynajian et al., PNAS 107, 10383 (2010). 

[3] J. Kondo, Prog. Theor. Phys. 32, 37 (1961). 

[4] A.C. Hewson, The Kondo Problem to Heavy Fermions 

(Cambridge University Press, 1993). 
[5] M. B. Maple et al., J. Low Temp. Phys. 99, 223 (1995); 

A. Schroder et al., Nature (London) 407, 351 (2000); 

J. Custers et al., Nature (London) 424, 524 (2003); P. 

Gegenwart, Q. Si, and F. Steglich, Nature Phys. 4, 186 

(2008); for reviews see: G. R. Stewart, Rev. Mod. Phys. 

73, 797 (2001); H. v. Lhneysen, Rev. Mod. Phys. 79, 

1015 (2007). 

[6] S. Doniach, Physica (Amsterdam) 91, 231 (1977). 

[7] Q. M. Si et al., Nature (London) 413, 804 (2001); P. Cole- 
man et al., J. Phys.Condens. Matter 13, R723 (2001); T. 
Senthil, S. Sachdev, and M. Vojta, Phys. Rev. Lett. 90, 
216403 (2003); I. Paul, C. Pepin, and M. R. Norman, 
Phys. Rev. Lett. 98, 026402 (2007); Y.-F..Yang et al., 



5 



Nature (London) 454, 524 (2008). 
[8] M. Maltseva, M. Dzero, and P. Coleman, Phys. Rev. Lett. 

103, 206402 (2009). 
[9] J. Figgins and D. K. Morr, Phys. Rev. Lett. 104, 187202 

(2010). 

[10] S. Hershfield, J. H. Davies, and J. W. Wilkins, Phys. Rev. 

Lett. 67, 3720 (1991). 
[11] O. Ujsaghy, J. Kroha, L. Szunyogh, and A. Zawadowski, 

Phys. Rev. Lett. 85, 2557 (2000). 
[12] M. Plihal and J. W. Gadzuk, Phys. Rev. B 63, 085404 

(2001). 

[13] J. Merino and O. Gunnarsson, Phys. Rev. B 69, 115404 
(2004). 



[14] V. Madhavan, W. Chen, T. Jamneala, M. F. Crommie, 
and Ned S. Wingreen, Phys. Rev. B 64, 165412 (2001). 

[15] K. Nagaoka, T. Jamneala, M. Grobis, and M. F. Crom- 
mie, Phys. Rev. Lett. 88, 077205 (2002). 

[16] Y. -F. Yang, Phys. Rev. B 79, 241107 (2009). 

[17] M. Fogelstrom et al., cond-mat/ 1004. 1882 (unpublished). 

[18] U. Fano, Phys. Rev. 124, 1866 (1961). 

[19] P. Cheng et aJ., cond-mat/1001.3220 (unpublished). 

[20] Y. Dubi and A. V. Balatsky, submitted. 

[21] E. Abrahams, and P. Wolfle, Phys. Rev. B 78, 104423 
(2008); P. Wolfle, and E. Abrahams, Phys. Rev. B 80, 
235112 (2009) 



